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Abstract 

An overview is given of Chiral Perturbation Theory and various appli- 
cations to eta decays. In particular, the main decay r/ ^ Svr is discussed 
at the one-loop level, and estimates of the higher order corrections are 
given. The importance of and higher order effects in double Dalitz 
and T] — > 7r'^77 decays is pointed out. In all these cases the need for new 
experimental results is stressed. 



-'^ Talks presented at the Workshop on Eta Physics, Angstrom laboratory, Uppsala, October 
25-27, 2001; pubhshed in the proceedings (Physica Scripta Topical Issues T99 (2002) 34). 
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1 Introduction 



Eta decays at high precision, both in measuring rare decays and in precisely 
determining kinematical distributions in the more common decays, have a high 
potential to teach us about various aspects of the strong interaction. In the 
view of the new experimental results to be expected from WASA at CELSIUS in 
Uppsala and KLOE at DAPHNE in Frascati and the ongoing work of the Crystal 
Ball, an overview of the present situation as regards some of the main eta decays 
is appropriate. The importance of these decays and measurements can be judged 
from the fact that they were the main focus of several theoretical talks presented 
at this workshop, in addition to the ones by the present authors there is also the 
contribution by Holstein|l|] and by Ametller0. 

In this review we concentrate on chiral symmetry aspects for the main strong 
decay 77 — Svr and how it can be used to extract information on quark mass 
ratios. We will show how precise measurements of the slope parameters will lead 
to more refined theoretical predictions for the decay rate as a function of the 
quark mass ratio 

md-rriu' 

A discussion of the decay 77 — > Svr in Chiral Perturbation Theory (CHPT) and the 
use of dispersion theory to improve on the one-loop prediction is the main focus 
of this short review. In addition, a short overview is given of Dalitz (77 '-jl^l^) 
and double Dalitz (77 decays, as well as of the process r] 7r°77. 

We indicate why these are important quantities to be measured. 

We first give a very short overview of chiral symmetry and its incorporation 
using effective field theory methods - known as Chiral Perturbation Theory - in 
Sect. ^. The decay 77 ^ Stt is then discussed at the purely perturbative level in 
CHPT at order in Sect. |] and the estimates of higher orders using dispersion 
relations in Sect. Section || is devoted to the Dalitz and double Dalitz decays, 
and the decay rj n^jj is the subject of section |^. Our conclusions are presented 
in section |^. 

2 Basics of CHPT 

We cannot calculate all things we are interested in in the strong interaction di- 
rectly from QCD. E.g. scattering amplitudes at low energies cannot be expressed 
as an expansion in the strong coupling constant because: 

• At low energies the coupling constant becomes large. 

• We deal with bound states which require always to go beyond a simple 
perturbation series. 
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So how can we calculate e.g. the cross-section of Trvr-scattering in QCD? 
Two fundamental options are possible: 

Lattice: This can be done using the methods developed by LiischerP via QCD 
calculations at finite volume. The energy spectrum at finite volume depends 
on the TTTT interaction. By measuring the lowest two-pion energy level, the 
TTvr S"— wave scattering lengths Oq can be obtained via 

= 2m. - + 0{L-') , (2) 

where / denotes isospin. There are first lattice results available which make 
use of this relation Q. 

Effective Theory of QCD: Using the methods of effective quantum field the- 
ory (EFT) we can also evaluate cross-sections in QCD at low energies. We 
replace the QCD Lagrangian by the effective Lagrangian of Chiral Pertur- 
bation Theory (CHPT). We replace the quarks and gluons by the degrees 
of freedom which are relevant at low energies: vr, K, rj, p, n,. . . . This is the 
method most of this report is concerned with. 

Chiral Perturbation Theory does reproduce for appropriately chosen coupling 
constants the 5'-matrix elements of QCD0] and it allows to make very sharp 
predictions in some cases. As an example, we quote the recent result on the 
scattering length combination ctj) — 02 lU > 

-ag = 0.265 ±0.004. (3) 

For earlier work at tree level and at one and two-loops accuracy, see [0, |], 3- 

We will now present the basics which underly the calculations of tttt tttt, 
rj TTTTTT, etc. in this approach. 



2.1 Chiral Symmetry 

Consider QCD with 3 flavours, u, d and s. Gluon interactions do not change 
the helicity of a quark, only mass terms couple the left and right handed helicity 
states. So in the limit where m„ = = = 0, the left-handed world cannot be 
turned into the right handed world. Now if all the masses are zero, there is also 
no way to distinguish the different flavours of quarks and we can continuously 
rotate one into the other. The QCD Lagrangian is thus invariant under 




Vj 



H 



d \ ; H = L,R, (4) 
V ^ Jh 
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separately for both helicities, left (L) and right (R). The matrices Vh are general 
unitary 3x3 matrices. So we have a global^ 



SU{3)l X SUi3)R (5) 

chiral symmetry of QCD . This symmetry has 16 conserved currents and has 
thus 16 conserved charges. The vector charges Qy = Q^j^ + also annihilate 
the vacuum, but the axial charges Q\ = —Q\ + produce a change on the 
vacuum: 

g{,|o) = o gi|o)^o. (6) 

The symmetry of the QCD Lagrangian is not the symmetry of the vacuum. This 
is known as spontaneous breakdown of chiral symmetry or spontaneous 
chiral symmetry breaking ((S)xSB). 

The 8 spontaneously broken axial symmetries require the existence of 8 Gold- 
stone bosons which must be massless and pseudoscalar. They are the pions, kaons 
and eta. The Goldstone bosons also have zero interactions at zero energy. The 
expansion of CHPT is based on these facts. 

Now, the fact that the quark masses are present and we can also have nonzero 
energies allows nonzero values for masses and interactions. E.g. the pion mass 
is an intricate mixture of the current quark mass, explicit x^B, and the quark 
vacuum expectation value, spontaneous x^BW^. 



ml = -{mu + md){Q\uu\Q)/F^ . (7) 
Similarly the pion scattering lengths do not vanish |^ 

aO = J^ = 0.159. (8) 

These predictions are not exact, they are the first term in a series expansion of 
momenta and the quark masses. 

2.2 A systematic expansion 

The amplitudes are expanded in powers of momenta and quark masses. The price 
paid for the generality is that order by order we get more terms, 

>CCHPT = C2 + C^ + C^ + ... , (9) 

where contains terms with n-derivatives or equivalent. In amplitudes the 
derivatives turn into momenta such that the expansion converges at low momenta. 
As said earlier, quark masses are counted as two powers of momentum, this is 



^The U{1) parts of the groups play no role here, but see the contributions by Bass|10 
Shore p), Michael pi, KroU and Feldmannpl for their effects. 
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because Eq. (|^ sets a quark mass to the pion mass squared. The structure of 
the Lagrangians is fixed by chiral symmetry but the couphng constants are 
not. 

In the strong and semileptonic mesonic sector, these Lagrangians are known 
and the parameters have conventional names: 



£2: Fo,Bo ^ 
Cq: Ci, . . . , C90 WT\ 



These couphng constants (LEC's) are free in chiral perturbation theory, but are 
in principle calculable from QCD0. 

The Lagrangians >C2,4,6,... allow to calculate S-matrix elements order by order 
in powers of momenta and quark masses. The main point is that order by order, 
only a limited number of terms can contribute, and the form of the contribution 
is completely fixed by chiral symmetry. So an amplitude can be expanded in 
powers 

A = A2 + A4 + Ae + --- , (10) 

with 

A2: tree graphs with vertices from £.2- The form of the vertices is fixed by chiral 
symmetry. 

A4: tree graphs with vertices from £2 and one vertex from £4 or loop graphs 
with vertices from £2- 



This procedure is known as Chiral Perturbation Theory (CHPT). 

In loop integrals one integrates over all momenta, from eV to the Planck mass 
and beyond, yet we always use the vertices from the effective Lagrangian con- 
structed to reproduce the amplitude at low energies. This is not a contradiction, 
the high energy behaviour of the loop integrals is irrelevant in the sense that 
this is the part that is absorbed by the coupling constants. In this way, different 
regularizations can be used to cut-off the integrals and the difference in different 
ways of doing it is absorbed by using different values of the coupling constants. 
The final result is always independent of the regularization used. In practice 
we nearly always use dimensional regularization since it preserves chiral sym- 
metry throughout the calculation reducing the need for technically complicated 
subtractions to restore symmetries. 

•^For some attempts to determine them from lattice QCD, see Ref. p8[ . 
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Ptt- 



Figure 1: The kinematical variables in the decay rj vr+vr vr". 
Many apphcations exist, see e.g. the miniproceedings of the Bad Honnef 



meetings and the review articles[lS] as well as the lectures EO]. 



3 ?7 ^ TTTTTT at order 

3.1 Kinematics and Isospin Relations 

The 1] is an isospin singlet and a pseudo-scalar. The decay into two pions is 
forbidden by CP. CP violation in rj decays is discussed by C. Jarlskog and 



E. Shabalin|21[] and J. Ng|^ in these proceedings. Three pions in an angu- 
lar momentum configuration cannot have isospin zero as well, but isospin 1 is 
allowed. This decay thus has to proceed via isospin breaking effects. Electro- 
magnetism is known to play a fairly minor role here as discussed below, except 
via the kinematical effects due to the charged and neutral pion mass difference. 
This decay thus goes primarily through the strong isospin breaking part of QCD 

= -^{ruu - ma) (uu - dd^ . (11) 

This itself has isospin 1 and there is thus to lowest order in isospin breaking a 
relation between t] —>■ Tr+Tr^Tr" and rj — >■ tt^tt^tt^. 

The kinematics is depicted in Fig. [l|. 

The variables used in the remainder are 



5 

U = {p^+ +P-„of = {Pr^-p^-f , (12) 

s + t + u = w?^ + 2ml+ + mlo = Ssq ■ (13) 



which satisfy 
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The last equality is the definition of sq. 

The decay amplitude itself can be written as 

(7r°7r+7r"out|?7) = i {27c'f 5^ (p^ — — p^- — p^o) A{s, t, u) . (14) 

Charge conjugation requires this amplitude to be symmetric under the inter- 
change of TT"*" and 7r~ so we have 

A{s,t,u) = A{s,u,t) . (15) 

Using isospin as discussed above and labeling the three momenta as pi,P2,P3 and 
Si = [Pr] — Pif , the amplitude for the neutral decay 

(7r°7r°7r°out|r/) = i (27r)' 5^ (p, -p^-p^- Ps) A{s,, ^2, s^) (16) 

satisfies 

^(Sl, S2, S3) = A{si, S2, S3) + A{S2, S3, Si) + A{S3, Si, S2) , (17) 



see e.g. [23] for a detailed derivation of this result. Experimentally, there are two 



main numbers to be kept in mind|^. The decay width from an overall fit 



T{r] ^ 71+ TV- 71^) = 271 ±25 eV , (18) 
with a scale factor of S* = 1.8, and the ratio of neutral to charged decays 

r = -^ = 1.404 ±0.034, 19 

with a scale factor of 5 = 1.3. The presence of these scale factors is an indication 
of mutually incompatible experiments. 

3.2 Lowest order: 

The lowest order CHPT contribution from the quark mass difference was derived 
in the sixties using current algebra methods and is 



As announced, the decay amplitude is proportional to m„ — itlii. We can now 
rewrite the amplitude in terms of Q, a combination of quark masses which is 
experimentally more easily accessible. 



2 "2 

mt — m 



Q = z:i— > (21) 



— mt 
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where rrii stands for the up, down or strange current quark mass, and 

rh = ^{rriu + rrid) . (22) 
The decay rate can then be written as 

. . ^^Ki 2 2 A M{s,t,u) 

^(^'^'-) = gi^(-^--x)^7^' (23) 



with at lowest order 



3s — 4m^ 



M{s,t,u) = — (24) 



as 



It can alternatively be rewritten in terms of 

R = — (25) 

nid - ruu 

A{s,t,u) = -^-±^{3s-4ml). (26) 

So if we have a good theoretical description of M(s, t, u), we get an accurate 
determination of Q, since 

r (r/ ^ Svr) oc |Ap oc Q-^ (27) 

Due to the high power, the error on Q can be made much smaller than its 
determination from meson mass ratios, 



because this relation has much larger electromagnetic corrections] 26 



In order to show the accuracy which can be obtained, we can input two 
different values of Q or R. We first use the value 

R = 40.8 ±3.2 [27], 
r(r7 ^ 7r+7r"7r°) ^ 66 eV, 

r = 1.51. (29) 

The second one is where we use instead the fit done at order for the quark 
masses from the meson masses and the estimate of the kaon electromagnetic mass 



difference including higher order effects [p^ , 

Q = 20.0 ±1.5 or i? ^ 29 , 
r(r/ ^ 7r+7r"7r°) ^ 140 eV. (30) 

It is clear that a good understanding of M(s, t, u) can lead to a very accurate 
determination of Q and R. 
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3.3 Electromagnetic effects 



The lowest order contribution does not contribute directly to the decay rate. 
It only comes in via the electromagnetic mass difference of pions. This is very 
relevant for the actual values of the kinematical variables s, t, u in terms of the 
measured momenta, but this can be included in a simple fashion. That this decay 
could not be explained by its electromagnetic contributions has been known for 
a long time already p9|. 

The calculation of the purely electromagnetic contribution to the decay has 
also been pushed to one- loop level |3^] and the total effect remains small. 



The decay with an extra photon in the final state has been estimated as well. 
This could have had potentially a large impact, since it can proceed without extra 



isospin breaking. Fortunately, both a simplified analysis using Low's theorem[^ 
and a full analysis at one loop in CHPT|^2| lead to a rather small result of the 
order of 1% of the decay width of r/ ^ Svr. 

The correct incorporation of the electromagnetic corrections also tells us which 
value for the pion mass should be introduced in the expression for the amplitude. 
It is the same for both decays and is the neutral pion mass but we need to include 
the correct physical masses in the phase-space integrals. This follows from adding 
to £2 the effective electromagnetic term 



with 



C2-^C2 + C{QUQU^) , (31) 

iQUQU^) = (tt+tt- + K^K-) + 1^^^+^- + • • • • (32) 

There is no direct contribution to 77 Svr. The first term leads to an equal 
mass shift for the charged pion and charged kaon in the chiral limit. This fact 



is known as Dashen's theorem ||33||. The second term generates a contribution 
through mixing. This may appear to be negligibly small, because it is of order 
e^{md — irtu)- However, this matters in the ratio r. Taking these corrections into 
account tells us that the pion mass appearing in Eq. ( ^61) is the neutral pion 
mass. 



3.4 Conclusions from the tree level 

• The total rate is off by a factor of about 4, depending on the precise inputs 
used. 

• The ratio r is near its experimental value. The reason for this is discussed 
below. 

• Second order isospin breaking effects are important in the phase-space. 
Since the latter is small, small terms are important there - they significantly 
lower r. 
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The effective Lagrangian method is a very efficient tool to perform these 
calculations. 



3.5 One-loop contribution 

The amplitude can be expanded as 



A{s, t, u) = A2{s, t, u) + Ai{s, t, u) + 



(33) 



A2 has been calculated in the previous subsections. To evaluate ^4(5, t,u), one 
needs to evaluate one-loop diagrams with C2 and tree diagrams with £4. This 
was done in Ref. The decay rate can now be written in the form of Eq. (^3]) . 
From the expression in Ref. [34|, it can be seen that the only free parameter 
coming in is L3. This parameter can be determined from 7r7r-scattering|T^ or 



more accurately from K^./^ decays ||2^, |35 



What one numerically finds is a very large enhancement over the lowest order 
expression, 

j dLIPS\A2 + Ai\ 



I 



dLIPSlAo 



2.4, 



(34) 



where LIPS stands for Lorentz invariant phase-space. This enhancement has as 
one major source the large S'-wave final state rescattering as was expected, see 
Ref. 



36 



and references therein. In Sect. ^ we will come back to estimates of 
this contribution to higher orders. This contributions is about half of the total 
enhancement found in Ref. 153 . 



3.6 r and slope parameters 

Can r be very different from 1.5? The lowest order amplitude has a zero for 
s = (4/3)m^. Higher order corrections are not expected to remove the zero and 
also not to move it very much. We can thus expect that the amplitude will remain 
roughly of the form 

A{s,t,u) ^b[s-Xml) , (35) 

with X about 4/3. Given this form we see from Fig. ^ that unless the zero moves 
very much, the ratio r will remain close to 1.5. The reason for this relation 
between X and the ratio r is the isospin relation between the amplitudes for the 
neutral and charged decays. 

At 0{p'^) the ratio r is changed from the tree level prediction. It is now 
lowered to 1.43, bringing it into nice agreement with the observed ratio. Again, 
the relation between X, or more generally the slope measurements in the Dalitz 



plot, and r follows from isospin. This is discussed in detail in Ref. |34 
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Figure 2: The dependence of r on X as defined in Eq. 
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Figure 3: Leutwyler's ellipse and the constraints from mass differences at 
Taken from Ref . . 



3.7 Q at order 

The quark mass ratio Q determines the major semi axis of Leutwyler's ellipse in 
the rris/md and rriu/md plane| 



mi 1 m 
— - + 



Q2 ml 



1. 



(36) 



The situation at order is shown in Fig. ^. The role that a better understand- 
ing of ?7 — s> 37r has in this context is shown in Fig. ^, where at order the 
measurement of — 37r is compared with the determination from the K^-K^ 
mass difference. The different points on the graph correspond to various esti- 
mates of the electromagnetic part of the mass difference. Note that estimates of 
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Figure 4: 
Ref. 1271. 



The importance of 77 — Svr for the measurement of Q. Taken from 



the contributions [^] indicate possibly even lower values of Q than are shown 
here. 



4 Tj ^ TTTTTT beyond order 

4.1 The method and equations to solve 

We can now try to estimate even higher orders using dispersion relation tech- 
niques. This should catch a large part of the higher order corrections, because 
the large vrvr rescattering part was already a large part of the one-loop prediction. 
Unfortunately, for a general three-body decay, this is not such an easy under- 
taking. There are several subtleties involved here, as compared to the simple 
case of a two-body final state, because one can have imaginary parts not only 
from the scattering, but also from the unstable particle. The latter difficulty 
can be avoided by studying instead the scattering process nr] — > vrvr with an rj 
mass such that the eta cannot decay. Then we only have tttt singularities in the 
various channels to deal with. This procedure is how one can restrict oneself to 
the influence of two-body scattering only. The general framework is known as 
Khuri-Treiman equations [^] and was used in Ref. [38| to study higher order 



corrections to eta decay. 

Here we will restrict ourselves to a simpler formalism that has been very 
successful in the study of tttt scattering[^. The underlying observation is that 
there are no imaginary parts from rescattering with angular momentum larger 
than or equal to two up to p^. This together with the fact that s + t + u is 
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constant shows that to order the amphtude ( p3[ ) can be written as 
M{s,t,u) = Mo(s) + (s-M)Mi(t) + (s-t)Mi(t) + M2(t) + M2(M)-^M2(s). (37) 

The functions Mj{s) correspond to isospin / rescattering in the two particles 
whose kinematics is described by s. The spht into these functions is not quite 
unique. There is some ambiguity in the distribution of the polynomial terms over 
the various Mj, because s + t + m is constant. The split (^) is useful because: 

• Dispersion problems are now reduced to a set of functions of only one 
variable much simplifying the analysis. 

• This split is fully correct to two-loop order. 

• The main two-body rescattering corrections happen inside the functions 
Mi. 



This method was used to estimate higher order corrections to rj Sir in Ref. 

A detailed description can be found in Ref. ||23|. With the approximations used 

in Ref. the two methods are in fact entirely equivalent. 

So we use here the form (p3|) and discuss the process tti] — > tttt with the final 
two pions in isospin I for an eta mass below the three pion threshold. This leads 
to the dispersion relations 

. N 1 /""^ , , ImM/(s') , , 

Mj{s) = - ds'- (38) 



up to subtractions, see below. To be more precise, the integrand contains the 
discontinuity 

ImM/(s') — ^ discM/(s) = ^ {Mi{s + ie) - Mi{s - ie)) . (39) 

We will need to subtract to be able to compare with the one-loop expression. 
Checking the high-energy behaviour of the one-loop expressions leads to three 
subtractions for Mq and M2, while two are sufficient for Mi. So the equations we 
will need to solve for Mj are 

s^ r ds' discMnfs') 



s f 

Mo(s) = ao + bos + Cqs'^ -\ / 

vr J 

s^ f ds' discMi(s' 
vr J 



s'^ s' — s — ie 



s'^ s' — s — ie 



/ N , s'^ f ds' discM2(s') 

M2{s) = a2 + b2S + C2s' + - I — -^—ii-i. (40) 



TT J s s — s — le 
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This leads to eight constants which should be determined. But the ambiguity 
in the precise choice of the M/ is visible here too. Up to the same order in the 
s, t, u expansion we have only four parameters 

M{s,t,u) = a + bs + cs"^ - d{s^ + tu) , (41) 

so there is some freedom in the choice. 

We can form some more convergent combinations at one-loop as well. Com- 
paring expressions one sees that 

M^{s) + ^M2{s) (42) 

and 

,M,M + M,W + ,^l^|_lIZ(5i|!) (43) 

F^im^^ - mi) 

need only two subtractions. This allows to determine c and d in terms of L3 and 
integrals over the discontinuities via 

4 1 /■ ds' ( 4 
c = Co + 3C2 = - y ^ |discMo(s') + -discM2(s') 

d = + - / ^{s'discMi(.0 + discM2(.')}. (44) 

So once we know a, b and L3, this forms a set of integral equations that can 
be solved. Now we have to put in the knowledge of the vrvr phases. In a single 
channel problem without singularities in the cross channel, this problem is simple 
and was solved by Omnes long ago. Here we have cross-channel singularities but 
they are again tttt scattering in an isospin 0,1,2 state. The discontinuities thus 



obey[41 



discMj(s) = sm6j{s)e~'^''^''> {Mj{s) + M/(s)} , (45) 
with Ml a consequence of the singularities in the t and u channel. They satisfy 

Mj(s) = J2fdcos9 cos" dcnirMp{t) , (46) 
where t is the appropriate t for a scattering of angle 6 with center of mass energy 



squared of s. The coefficients Cnir can be found in ||23|, ^|. 

So at least we have a principle solution. There are several technical complica- 
tions involved. After analytically continuing m,^ to its physical value, extra sin- 
gularities appear, and we need M{s) for values of s outside the physical domain, 
so singularities in the relation of t with s and cos 9 appear. The integration path 
needs to be chosen carefully to avoid these extra singularities. Having solved this 
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problem, one finds that the solution to the equations is not unique. The problem 
is that homogeneous equations of the type 



1 /■ (^s' 

Mi(s) = - / sm6i(s)e-'^''^'^ {mAs) + MAs)} (47) 

IT J s — s — ^ J 

have nontrivial solutions - we need to pick out the physically correct one. The 
solution is to write a dispersion relation for a related function which only 
has one solution. We remove the singularity in the direct channel by dispersing 
instead 

-,W = ^. (48) 

with 

^ , , ( s r ds' Sj(s') ] , , 

nj{s) = exp - / — — . (49) 



TT J s s — s — te 



Similar arguments to above lead then to the dispersion relations with a total of 
4 subtraction constants for the m/. They are, rewritten in the M/, 





is) 




» 


Mil 


is) 


fill 


» 




is) 


n^i 





sm6o{s') Mo{s') 



"0 + /^os + 7oS^ + — / ds' 

IT J \ilo{s')\s''^{s' — s — le 

sin 6 i{s') Mi{s') 



n , s f , smdi[s 



— s — is)'' 



^ r sin 6,{s')M,{s') 

" TT J \n2is')\s'^is'-s-te)- ^ ' 

So basically we now have to find a parametrization of 5o.i,2(s), the vrvr phase- 
shifts, and determine the subtraction constants. After that the above equations 
can be solved numerically. Two of the constants can be determined similarly to 
c and d above, leading to 

70 ~ 0, 

ft . Jjizll^. ,51) 

The values of these two parameters compared with the one-loop expressions are 
basically independent of the matching points. The other two are discussed in the 
next section. 



4.2 Results 

The two papers |^ and basically only differ in the choice of the subtraction 
procedure. As explained above, the formalism is quite different but the underlying 
approximations and assumptions are the same. 
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Figure 5: The result of Ref. ||4^ for ReM (s, t, u) along the line t = u. 



Reference |^0| looks at the lowest order expression for 



M{s,t,u) = '-^^^, (52) 
— 

and notices that we can determine ao and /Sq from the place of the zero, the Adler 
zero happening at sa, and the value of the slope in s at s = sa- 

Reference instead uses several values of s, t, u to fit the dispersive approx- 
imation to the one- loop result. The resulting divergence in results is taken as an 
estimate of the theoretical error. The final result for the enhancement over the 
one-loop result is quite similar in both papers, but a closer look at the intermedi- 
ate results shows some discrepancies which need to be understood. The one thing 
which can be compared are the published plots of ReM(s, t, u) as a function of s 
along the line t = u. These are shown in Fig. |^ and Notice that while both 
have a significant enhancement, the slope is quite different. The results of |5D| are 
always above the one-loop result while the results of [^] cross the one-loop result 
at the end of the physical region. Given the similarity of the methods and the 
inputs used, this type of discrepancies needs to be understood, but it is obvious 
that a better experimental determination of the slopes can distinguish between 
both calculations. It is therefore very important that all slope parameters, see 
below for the definitions, are well measured. 

As a simple estimate of the uncertainties involved, we have used a simplified 
version of the equations (0) where we set the functions M/ = 0. So we only 
consider rescattering effects in the direct channel. We then fix the values of ao 
and Pq, considered real, to the tree level slope and Adler zero. In Fig. |^ the tree 
level is shown (tree) and the dispersive improved tree level (dispersive tree). The 
two dispersive improvements of the one-loop result shown are using the position 
of the one-loop Adler zero, and either the absolute value of the slope (dispersive 
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Figure 7: The value of ReM(s, t, u) in a simplified version of the dispersion 
relations with several different estimates of the subtractions. and (3q are de- 
termined from the tree level (dispersive tree) the position of the one-loop Adler 
zero, and either the absolute value of the slope (dispersive absolute) or its real 
part (dispersive real) at one-loop. The tree level (tree) result is shown for com- 
parison. 



absolute) or its real part (dispersive real) at one- loop to determine and [3q. It 
can be seen that these predictions lead to different slopes in the physical regions, 
thus allowing the theoretical subtraction procedure to be experimentally tested. 
Notice that we have shown only one variation in the phase-space. The variation 
over all of phase-space can of course be used similarly. 

The total enhancement found is fairly small, about 14% in Ref. for the 
total rate and 5% in the center of the Dalitz plot in Ref. |40], with a total rate 
enhancement in very good agreement with Ref. [^. 



4.3 Dalitz plot expansions or slope parameters 

Once the amplitude is known, one can easily expand the Dalitz plot distributions 
around the center of the Dalitz plot. Notice that some care is needed to precisely 
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Figure 8: The Dalitz plot or phase-space of rj tt+tt^tt'^ decays. Shown are the 
boundaries of the physical region, the variation with x and y, see text, and the 
point s = t = u. 

define the center of the Dalitz plot, given the second order isospin breakings 
discussed earlier. The point s = t = u = Sq does not quite coincide with the point 
where all kinetic energies are equal for the charged decay. The usual definition 
of the slope parameters is in terms of the kinematical variables x and y, 

= rrir, - 2m+ - m° . (53) 

Tj stands for the kinetic energy of the pion of charge i. In Fig. ^ we have shown 
how s and u depend on x and y. The crosses correspond to x = and y equal 
to —1, — 0.9,. . . ,1. The X correspond to y = and x = —1,0.9,...,!. The 
filled square indicates the place where s = t = u. The difference with the point 
X = y = Q is again an indication of the isospin breaking effects in the masses 
discussed earlier. 

The measured distribution is now parametrized as follows. For the charged 

decay, one takes 

l + ay + hy'^ + cx^ , (54) 
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a h c g 


tree 
one-loop p4[ 
dispersive |38| 
tree dispersive 
absolute dispersive 


-1.00 0.25 0.00 0.000 
-1.33 0.42 0.08 0.03'^ 

-1.16 0.26 0.10 -0.014 0.028 

-1.10 0.31 0.001 -0.013 
-1.21 0.33 0.04 -0.014 



Table 1: Theoretical results for the slope parameters of the various approxima- 
tions discussed in the text. The spread gives an indication of the uncertainty, for 
g it is about a factor two, even if the central values agree much better, see e.g. 
the discussion in |^8|]. "Ref. does not quote a numerical value for g. We have 



obtained g = 0.03 from the representation of the amplitude provided in pl|. 



normalized to the center x = y = oi the Dalitz plot, and 

l + g{x^ + y') (55) 

for the neutral decay. It is hard to compare with the experimental results in 
the review of particle properties [Q. They simply state that the assumptions 



made in the different measurements are not compatible. The problem is that 
making assumptions on the values of the quadratic slopes significantly alters the 
fitted result of the others. In table |l| we show the theoretical results of tree level, 
one-loop, the dispersive results from |^ and the results of the tree dispersive 
and absolute dispersive simplified estimates discussed above. Note that the one- 
loop predictions of the neutral slope g is not zero, contrary to the statement 



made in Refs. ||42| , ^]. The loops themselves do give a contribution, whereas the 
contributions from the tree and tadpole diagrams at happen to vanish. The 

precise value is dependent on the procedure used and is sensitive to fairly small 
changes. The number quoted corresponds to the precise way the expressions are 



given in We plan to come back to this issue in a later publication. 

Measurements as mentioned earlier are difficult to compare, but several new 
results have been obtained recently. In particular, the new Crystal Ball value for 
g was announced at this meeting [|^. Note that the experiments usually quote 
a = g/2. For completeness the known experimental results are given in table ^ 
for the charged decay and table |^ for the neutral decay. 

Notice that while the overall agreement is all right, there is a problem in 
obtaining a value of g which is compatible with the newer experiments. The 
theoretical evaluation of g suffers from large cancellations. 
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a h c 


Layter |43| 
Gormleyp^ 
Crystal Barrel 45 1 
Crystal Barrel 


-1.08 ±0.14 0.034 ±0.027 0.046 ± 0.031 
-1.17 ±0.02 0.21 ±0.03 0.06 ±0.04 
-0.94 ±0.15 0.11 ±0.27 
-1.22 ±0.07 0.22 ±0.11 0.06 fixed 



Table 2: Experimental results for the slope parameters in the charged decay. 
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Alde|7 
Crystal Barrel 48 
Crystal Ball IS, 42 


-0.044 ±0.046 
-0.104 ±0.039 
-0.062 ±0.008 



Table 3: Experimental results for the slope parameter in the neutral decay. We 
have added in quadrature the statistical and systematical errors quoted in . 



5 Dalitz and double Dalitz decays 



This section is a shortened version of the work presented in Ref . [£0| . The under- 
lying question is how well do we understand the couplings of the r] to two photons. 
In table ^ we show the (expected) branching ratios for all of the decays with two 
photons. The decays with one lepton-antilepton pair are known as Dalitz decays, 
those with two pairs as double Dahtz decays. These decays allow us to study the 
off-shell structure of the r/77 vertex as shown in Fig. ^. The question is: how 
does this form-factor F{q\^ g|) behave over the entire range of off-shell masses 
for the photons? Is it of order 1, w.r.t. to on-shell photons, can it be described 
by adding vector meson propagators in the photon legs? One would also like to 



Decay 


Branching ratio 




39.33 ± 0.25 % 


Tj — > 76+6" 


(4.9 ± 1.1) ■ 10-3 




(3.1 ±0.4) ■ 10-^ 


7] —>■ e+e^e+e^ 


~ 6 • 10"^ 


T] e^e^ jj,^ jJL^ 


~ 2 • 10-6 


1] fi^^^fi^fi^ 


~ 10-8 



Table 4: The measured or expected branching ratios for various decays of rj to 
two (possibly off-shell) photons. 
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Figure 9: The 7777 vertex with photons of mass squared and 
know the approach to the perturbative QCD result (see e.g. Ref. |l5l[| ) 

F{-Q',-Q') = ^^ + ... ,Q'-.oo. (56) 

This type of form-factors are one of the places where experiments are possible 
for the same underlying physical process over the entire regime, ranging from the 
fully perturbative to the fully nonperturbative one. The large and intermediate 
values of and Ig^l can be studied in tagged photon-photon coUisions while the 
small timelike values can be studied in r^-decays. 



To study how well r^-decays can see variations in the form-factor, in Ref. |^ 
four trial form-factors were investigated, 

Fiqlq'2) = 1, (57) 

4 

(58) 



[m\ - (H) {ml - ql) ' 



mir , , 

(59) 



{ml - ql - ql) ' 
< - (4vr^F|/3) {ql + ql) 
{ml - ql) {ml - ql) 

The first form is included to check how well a deviation from a pointlike eta 
can actually be measured. The second one is the standard double vector meson 
dominance model, while the third was a variation that reproduced single vector 
meson dominance, while numerically approaching Eq. (|56D quite well. The last 
form was suggested in Ref. , see also Ref. |^ , Knecht and Nyffeler in Ref. |^ 



as well as the contribution by Feldmann and Kroll to this conference |]T3[. One 
important reason why we like to know this form-factor better is that pseudoscalar 
exchange is the major contribution to the hadronic light-by-light scattering part of 
the muon anomalous magnetic moment as studied in Ref. |Q. The contribution 



of the exchange of vr'', and 77' using the different form-factors to = {g^ — 2)/2 
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Decay 




053) (60) 




2me 
50 MeV 
200 MeV 
300 MeV 
400 MeV 


1.62-10-^ 1.67-10-^ 1.66-10-^ 
4.6-10-3 5.1-10-3 5.0-10-3 
8.6-10-4 11.5-10-4 10.9-10-4 
2.2-10-4 3.6-10-4 3.3-10-4 
3.0-10-5 6.3-10-5 5.6-10-4 






5.5-10-4 7.7-10-4 7.3-10-4 



Table 5: The branching ratios for the single Dalitz decays relative to — 77 as 
a function of the lower limit of the lepton pair mass for the different form-factors 
defined in the text. 



Decay 


^^e+e Iniin 


m m m 


T] — > e'^e^fi^fi^ 


2me 
50 MeV 
200 MeV 
300 MeV 


5.6-10-6 5.6-10-6 5.3-10-6 
11.6-10-^ 11.7-10-^ 10.9-10-^ 
5.0-10-^ 5.2-10-^ 4.7-10-^ 
2.8 - 10-1° 2.9 - 10-1° 2.6 - 10-1° 



Table 6: The branching ratios for the double Dalitz decay rj e^e^fi^jj,^ 
relative to ^ 77 as a function of the lower limit of the electron positron pair 
mass for the different form-factors defined in the text. 



is g 7.9 - 10-1° for form-factor (||), 9.7 - 10-^° for (H) and 10.9 - 10-^° 
for (|60D . It is thus clearly important for future measurements of that this 
form-factor becomes experimentally more constrained. 

In tables ^ to |^ we show the results |^ for the various decays as a function 
of the electron-positron mass lower bound as well the muonic cases for all of 
phase-space. It turns out that it is surprisingly difficult to see the difference 
between the various models of the form-factors in 77 decays, but some constraints 
are possible. The form-factors which are different in single Dalitz decays should 
be easily distinguishable, but more subtle differences, as the one between ( ^8| ) and 
(|59D - which only become visible in double Dalitz decays - will be difficult to see. 
Given the importance of these results, the measurements should be performed. 
The differences between the effects of the the form-factors in (0) to ( |60D is 0{p^). 
Form-factors (|58D (|59|) only differ at 0{p^). The differences are a good indication 
of the importance of the various orders. 



6 rj ^ TT 77 

This decay is only touched upon in this report. A more comprehensive discussion 
can be found in the review by Ametller[^ and in the references. The main theory 
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Decay 


^^e+e mill 


m m m 


T] — >• e^e^e^e- 


2me 
50 MeV 
100 MeV 
200 MeV 


6.7-10-^ 6.7-10-^ 6.6-10-^ 
4.4 • 10-6 4.4 ■ 10-6 4.3 ■ 10-6 
7.4 • 10-^ 7.5 ■ 10-^ 7.2 ■ 10"^ 
2.7-10-9 2.8-10-9 2.6-10-9 






9.2-10-9 9.4-10-9 8.5-10-9 



Table 7: The branching ratios for the double Dalitz decays rj jj,'^ jj,' jl^ jj,^ and 
7] — > e^e~e^e~ relative to ?7 ^ 77 as a function of the lower limit of the electron 
positron pair mass for the different form-factors defined in the text. The last line 
was not published earlier. 



paper underlying this decay is and the main interest is that this decay is a 
window on rather high order corrections in CHPT. The O(p^) result is small for 
several reasons. The amplitude has two possible Lorentz structures. In the gauge 
ej - Qj = 0, with ej and qi the polarization vector and momentum of photon i, the 
amplitude can be written as 

M(?7 7r°77) = Aqi - g2 ei - 62 - 5 (ei - €2 P,, - qi Pr, ■ q2 + P,j ■ ei Pr, - €2 gi - ^2) • 

(61) 

At 0{p'^) only the A part of the amplitude is nonzero and it is small. The pion 
loop is suppressed by isospin breaking and the kaon loop is suppressed by 4m|. 
and an extra factor of 1/12 in the integral. The contributions are 

vr-loops 0.84 - 10-^ eV , 
/sT-loops 2.45 - 10-^ eV , 

sum 3.9-10-^eV, (62) 



rather far from the old experimental value of 0.84 ± 0.18 eVp^. Newer exper 



imental limits come from Novosibirsk ||56|| and the Crystal Ball[^. The latter 
gives 0.42 ± 0.14 eV. The result of Eq. ( |62D is therefore more than two orders of 
magnitudes below the measured width. 

We can also add vector meson exchange as depicted in Fig. This starts 
at (9(p6) and has contributions to all orders. Restricting the exchange of uj and 
p to p6 terms only leads to r^6^ = 0.18 eV and if all orders are kept 

rVMD = 0-31eV, (63) 

much larger than the suppressed contributions. In fact, there is one more 
contribution with well determined signs which appears at 0{p^). Here there is 
no isospin suppression of the pion-loops any longer nor any factors like the 1/12 
at O(p^). As a result the p^ contribution from the double WZW vertex one- 
loop diagram is as large as the p'^ one. Given the structure of the amplitude, 
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Figure 10: The vector meson exchange contribution to — > 7r°77. 

interference effects are quite strong so that if we put in the loops plus the VMD 
estimates we obtain the estimate 

r^0.42eV. (64) 

We can now add the contribution from and 02 exchange, the absolute values 
are experimentally known but the signs are free. This leads to 

r ^ 0.42 ±0.20 eV. (65) 

in very nice agreement with the preliminary Crystal Ball result. Studies of distri- 
butions in the decays will allow to disentangle various contributions to this decay. 
It should be noted that studies of this decay mode have also been performed in 



the ENJL modelfS^] and by adding other resonances than the ones referred to 



above[BSl. There also exist some recent quark model studies! 



7 Conclusions 

We have discussed several eta decays from the point of view of Chiral Perturbation 
Theory and various possible enhancements. The main conclusion is that more 
and precise measurements are very much needed. 

7] ^ : This decay has the potential to deliver accurate quark mass ratios. 
For this we need better measurements of the slope parameters in the Dalitz plot 
and a better determination of the ratio r between the neutral and charged decay 
mode. These will allow to put the theoretical calculations of the decay rate on 
the experimentally verified footing needed to obtain accurate quark mass ratios. 
A possible discrepancy with the slope g in the neutral decay should be checked 
theoretically and further experimental confirmation would be welcome. 

rj — *• 7*7* : The double off-shell form-factors are needed in these double Dalitz 
decays and provide a useful insight into short-distance long-distance transitions in 
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QCD. They are also needed input in precision calculations for the muon anoma- 
lous magnetic moment. It turned out to be surprisingly difficult to see these 
differences in ?7-decays, but it is an important measurement to be performed. 

rj —>■ n^'-f'-f : The existing discrepancy with theory seems to be on the way to 
be resolved but we need con&mation of this and the study of distributions will 
allow different contributions to this decay mode to be distinguished. It is a decay 
mode where the usually dominant contributions CHPT at tree level and one-loop 
are very small. It thus provides a rare window on higher order contributions. 
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